Localization in ruthenates: magnetic and electronic properties of Ca2_2^^Sr2:Ru04 



O 

o 
o 

(N 
> 

o 

00 



I 

o 
o 



> 
o 



o 
o 



o 
o 



X 



V.I. Anisiniov\ LA. Nekrasov\ D.E. Kondakov\ T.M. Rice^ and M. Sigrist^ 
^ Institute of Metal Physics, Russian Academy of Sciences- Ural Division, 620219 Yekaterinburg GSP-170, Russia 
Theoretische Physik , ETH-Honggerberg CH-8093 Ziirich, Switzerland 
^ Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

(February 1, 2008) 



The electronic structures of the metalhc and insulating 
phases of Ca2-iSra;Ru04 (0 < a; < 2) are calculated us- 
ing LDA, LDA+U and Dynamical Mean-Field Approxima- 
tion methods. For x = Q the ground state is an orbitally non- 
degenerate antiferromagnetic insulator. For 0.2 < a; < 0.5 
we propose a state with partial orbital and spin ordering. For 
X > 0.5 the observed Curie- Weiss paramagnetic metallic state 
possessing a local moment with the unexpected spin 5=1/2 
is explained by the localization of only two of the three Ru- 
4d-orbitals. 
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The isoelectronic alloy series Ca2-a:Sr3;Ru04 is a rare 
example of a metal-insulator transition with a multi-band 
electronic structure. Sr2Ru04 has been well studied of 
late because of its unconventional superconductivity and 
strongly correlated Fermi-liquid behavior in the normal 
phase The other end member, Ca2Ru04 is a local- 

ized Mott insulator, antiferromagnetically (AF) ordered 
and with a distorted crystal structure [^|8| . It is very in- 
teresting to examine how this isoelectronic system evolves 
between these two contrasting end members. 

Very recently Nakatsuji and Maeno (NM) reported 
transport, magnetic and crystallographic measurements 
on the intermediate phases iQ,^. Especially notewor- 
thy is the pronounced increase in paramagnetism in a 
metallic phase as x is decreased from pure Sr2Ru04 
{x = 2) down to a critical value Xc — 0.5. In this range 
2 > X > 0.5 (called Region III by NM), the susceptibility 
evolves to show free Curie behavior (i.e. x « C /T) at 
X = Xc with a Curie constant C corresponding to a spin 
5 = 1/2 per Ru'*+-ion. This result presents two puzzles: 
the apparent absence of exchange interactions in this con- 
centrated spin system and the unexpected reduction in 
the value of the local moment from the usual S = 1, 
for localized Ru^+-ions [Q . Further reduction of a; < a;c 
leads to Region II characterized by a magnetic metallic 
phase in which both a tilting of the RuOg octahedra and 
AF correlations between the S* = 1/2 moments appear 
in a metallic phase. Finally in Region I (0 < a; < 0.2) 
a metal-insulator transition occurs as the temperature 
is lowered. The end member Ca2Ru04 is an AF Mott 
insulator with the expected 5* = 1 local moment and a 
more strongly distorted crystal structure. This complex 
evolution of the electronic structure from a well-formed 
Landau-Fermi liquid in Sr2Ru04 to the AF Mott insu- 



lator in Ca2Ru04 must be interpreted in terms of the 
multiband electronic structure of these materials. 

In this Letter we report calculations of the electronic 
structure of Ca2-a;Sra;Ru04 using local density (LDA), 
local density with onsite correlations (LDA-I-U) and Dy- 
namic Mean-Field (DMFT) approximation schemes, in 
each range for compositions for which crystal structures 
are determined. In this way we track the evolution in the 
electronic structure from the metal at a; = 2 to the AF 
insulator at x = 0. 

Sr2Ru04 has the single-layered K2NiF4-structure 
(space group lA/mmm) with layers of RuOe- 

octahedra, slightly elongated along the c-axis and tetrag- 
onal local symmetry for a Ru'^^-ion. The band structure 
has completely filled 0-2p-bands, 4 electrons in Ru t2g 
4d-band and an empty e^-band. Due to the tetragonal 
symmetry, the xy-orbitals of the t2g-band are not equiv- 
alent to the (a;z, ?/2;)-orbitals. The crystal field splitting 
between them is small but the xy-orbitals 7r-hybridize 
with the 2p-orbitals of all 4 in-plane 0-neighbors while 
xz (j/z)-orbitals 7r-hybridize only with the 2 0-neighbors 
on the X-axis (y-axis). As a result the xy-band has 
two-dimensional dispersion and a bandwidth nearly twice 
that of the (xz, j/z)-orbitals (see Fig. 1). The measured 
effective masses, however, are typically enhanced by fac- 
tors of 3 to 4, indicating strong correlations. Note the 
volumes contained within the Fermi surface sheets show 
an approximately equal occupation of all 3 orbitals — a 
state we denote as (8/3,4/3) where {n(^a,i3),nj) are the 
occupation of the (yz, zx)- and xy-orbitals, respectively. 

Turning to the other end member, we note that the 
ionic radius of Ca^+ is smaller than Sr^+, so that a stabi- 
lization of the metallic state might be expected. However 
with the substitution of Sr by Ca in Ca2-a;Sr2;Ru04, the 
average (Sr,Ca)-ion size becomes too small for its space 
in the lattice of RuOg-octahedra. The lattice reacts by 
rotating and tilting the octahedra so as to decrease this 
space. Also the smaller size of the Ca^^-ion decreases 
the interlayer distance (c-axis lattice constant) which 
results in a smaller elongation of the RuOe-octahedra 
(for Ca2Ru04 the octahedra are even compressed 0). 
These structural changes lead to a smaller bandwidth 
due to the reduction of the Ru-O-Ru-bond angle from 
the ideal value of 180° and the compression of the octa- 
hedra changes the sign of the xy-{xz, yz) energy splitting. 

We consider the evolution of the electronic structure 
starting from the good metal, Sr2Ru04. Here the frac- 
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tional occupation of all 3 bands that cross the Fermi 
surface clearly inhibits Mott localization, since none of 
the Fermi sheets allows elastic umklapp electron-electron 
scattering at the Fermi energy. Yet the evolution from 
a; = 2 to a; = 0.5 clearly points to the partial local- 
ization of the 4d-electrons. Metallic strongly correlated 
electron systems have been successfully described by the 
recently developed DMFT |ll|,|l|l. We have used LDA 
to determine the band parameters of the non-interacting 
Hamiltonian and then treated the interactions within the 
DMFT. The effective Anderson impurity model in the 
DMFT scheme has been solved by the non-crossing ap- 
proximation (NCA) lH^. This leads to the combined 
method called LDA+DMFT+NCA @. 

The most puzzling aspect is the effective S ^ 1/2 
moment observed in the Curie- Weiss-like susceptibility, 
since the ionic value for the low-spin state of a Ru'*^-ion 
is S" = 1. In the absence of full crystallographic data for 
this region, we applied the LDA-I-DMFT-I-NCA scheme 
to the Sr2Ru04 structure, with U values varying from 1.0 
eV to 2.5 eV. The partially filled Ru-4d-bands derived 
from the xy and (yz, zx) orbitals do not hybridize with 
each other by symmetry. Since the effective bandwidth 
of the xy-band is considerably larger (Fig. la), a larger 
value of U is needed to drive the metal-insulator tran- 
sition for the xy-hand than for the (yz, za;)-band. The 
LDA-hDMFT-hNCA calculations showed that Uc ^ 1.5 
eV is a critical value for localization of the {yz, zx)- 
electrons while a larger value oiUc = 2.5 eV is needed for 
the xy-hand. These results imply that for intermediate 
values oiU (1.5eV<C/<2.5 eV), there will be a regime 
of partial localization of the 4d-electrons. In practice the 
RuOe-octahedra rotate when a; is < 1.5, reducing the n.n. 
hopping and narrowing the band. For constant U, Mott 
localization will appear first in the {xz,yz)-hand. This 
leads us to propose that the region x ^ Xc has a novel 
electronic structure illustrated in Fig. lb. Of the 4 elec- 
trons in the i2g-bands, 3 electrons (or 1 hole) are in the 
Mott localized (xz, yz)-bands and are the origin of the 
S — 1/2 local moments. The remaining 1 electron forms 
a half-filled but still metallic xy-hand. Integer occupancy 
of the different sub-bands is the key to achieving Mott lo- 
cahzation. Our LDA DFMT NCA calculations give 
a self-consistent solution with exactly (3, 1) occupancy. 

We now turn to the magnetic properties of the local 
5 = 1/2 moments. While a Kondo-type of interaction 
between the two bands (xz,yz) and xy can be excluded 
due to the absence of hybridization for symmetry rea- 
sons, an AF-type of RKKY interaction between the lo- 
calized spins is still induced via Hund's rule coupling. At 
the same time, however, Hund's rule coupling causes the 
ferromagnetic correlation through the double exchange 
mechanism. Thus, the two types of spin interactions me- 
diated by the itinerant electrons of the a;j/-band are likely 
to compensate each other, such that the exchange cou- 
pling between neighboring localized spins occurs mainly 
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FIG. 1. Density of states (DOS) for Sr2Ru04: a) The 
partial DOS obtained from the LDA calculation, b) The 
results of LDA-I-DMFT-I-NCA calculations for paramag- 
netic metallic state. The parameters for non interact- 
ing Hamiltonian were calculated using Sr2Ru04 crystal 
structure and Coulomb interaction parameter U=2 eV. 
The solid line is the DOS for (xz, y2;)-orbitals and dashed 
line for xy-orbital. (n{yz,zx).,'n,7:v) indicates the electron 
distributions. 



through superexchange processes in the (yz, zx)-band. 
The highly anisotropic hopping matrix elements between 
these orbitals, however, leads to an essential dependence 
of the superexchange interaction on the orbital config- 
uration of the minority spin electron (or single hole) of 
each Ru'*+-ion in the degenerate [xz, yz) bands. To gain 
insight into the possible form of orbital ordering we per- 
formed LDA+U calculations for the critical concentra- 
tion Xc = 0.5 ll|]. 

Cai.5Sro.5Ru04 has the space group l^xjacd |l^] 
which has no tilting of RuOg-octahedra, only rotations 
around the c-axis. The RuOe-octahedra are also elon- 
gated causing the a;y-orbital to be higher in energy than 
[xz, j/z)-orbitals leading to orbital degeneracy in a (3, 1) 
state. The LDA-I-U calculation gave an orbital order of 
an antiferro-type. The minority-spin electrons occupy al- 
ternating xz- and j/z-orbitals ((t^) = ±1) with a slight 
tilting of the orbital planes away from the c-axis. This or- 
bital order favors ferromagnetic (FM) spin alignment and 
this could explain the experimentally observed enhance- 
ment of susceptibility with increasing of Sr concentration 
toward the x = 0.5 

Specifically the localized orbital degrees of freedom 
are represented by the isospin r where the yz- and zx- 
orbitals on each site correspond to the states |-|-) and |— ), 
respectively, with t^|±) = ±i|±). We parameterize the 
n.n. superexchange interaction as 

= ^ Ja(Tj,-rj+a){Si ' Si+a + (1) 

with Jaiji, Ti+a) = Jili,a ' Ii+a,a- Wc define the orbital 
(2-component) vector operator Ii_a = {jl -I- gaT^ ,Tf -|- 
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ga'rf) where a = (1,0), (0,1) is the vector connecting 
n.n. sites and g = — a-y (r° = ^12x2)- This Hamil- 
tonian is reduced to the essential parts from the most 
general possible form. The parameters a, a', b and Ji 
are chosen to be consistent with the results of the LDA 
calculation. In particular, < 1 — a <C 1 arising from 
the reduction of the overlap for orbitals with their planes 
_L a. This leads to AF coupling, if the orbital states 
on a bond are the same (strongly reduced, if the bond 
is perpendicular to the plane), and to FM coupling for 
the configuration |+) (g) |— ). This competition between 
AF and F coupling suggests that there may be a strong 
cancellation due to fluctuations of the orbital occupation. 

At lower values of x we enter Region II (0.5 > x > 
0.2) characterized by a tilting plus rotation of RuOe- 
octahedra. Cai.8Sro.2Ru04 has a low-symmetry crys- 
tal structure with the space group P2i/c [|o), which 
can be obtained from the tetragonal lA/mmm structure 
by rotating and tilting of the RuOe-octahedra similar 
to pure Ca2Ru04 but with a smaller tilting angle 
There are now two types of in-plane oxygen ions and 
two types inequivalent of RuOg-octahedra. The ratio of 
the Ru-0 bond length for the apical and planar oxy- 
gens is larger than 1, so that the xy-level lies higher 
than the {xz,yz)-\evels leading to an occupation (3, 1). 
The single minority-spin electron is orbitally degenerate 
in {yz, zx)-band. We performed LDA+U calculation for 
this tilted structure and obtained a rather complicated 
orbital order. The ground state is an AF insulator. The 
minority-spin electron (1 per Ru-atom) occupy the or- 
bitals whose planes are in average directed along the a- 
axis (in tetragonal notation (110) direction). However on 
every one of the 4 Ru-atoms in the unit cell those planes 
are rotated from the a-axis by -1-20° ((r) = (it,0, u)) 
and +15° ((r) = (m',0, w')) on one layer and by -20° 
((r) = (-u,0,w)) and -15° ((t) = {~u',0,v')) on the 
next layer {u = 0.643, v = 0.766, u' ~ V3/2 and 
v' = 1/2). Also on one of the 2 Ru-atoms in every layer 
there is an additional tilting of the orbital plane from 
the long c-axis on 34° which is not represented within 
our model. The calculation of the easy axis using the 
second order perturbation theory for spin-orbit coupling 
gave the direction of magnetic moment as along the a-axis 
(tetragonal [110] direction) with a 28° tilt from the layer 
plane. The measurements of the uniform magnetic sus- 
ceptibility show a peak in the temperature dependence 
which is the most pronounced for the [110] direction of 
the magnetic field, in agreement with our LDA+U results 

i- 

To include the aspect of spin anisotropy we extend the 
interaction adding the following terms, 



T^'^ay VV-'i.a-'j-ha,a ^ -'i,a-'i+a,a *i+a ^ *i "t+a) 

(2) 

+ -^5Sf<+a] 

where further phenomenological parameters appear. Our 
previous discussion based on LDA indicated that crystal 
field effects yield the dominant bias for the orbital orien- 
tation. This is taken into account by an additional term 
Hcf = - Y.^,^,Q^'rt where is a local "Jahn-Teller" 
mean field which depends on the distortion introduced 
by Ca-doping and may show a complicated staggering as 
suggested by our LDA calculation. 

We now can connect the experimental results for x = 
0.2 with our model. Taking for simplicity Qi (x (1,0,0) 
uniformly plane of all orbitals is parallel to [110] and we 
obtain an effective AF anisotropic n.n. spin exchange of 
the form ^ ■ ^ ^ Jt^SiS^^^. Here the coordinates ^ are the 
principal spin axis in the exchange (1) and (2). Two axis 
(/i = x' and z) lie in the plane of the orbitals and the 
third axis (/i = y') lies orthogonal. Assuming J3 > 
and J5 K, J3 leads to an "XXZ-type" Heisenberg model 
where the Ising part (spin quantization axis along [110]) 
is weaker. Consequently, the uniform susceptibility is 
larger along the [110]-direction. The temperature depen- 
dence of the susceptibility for all directions show peak 
features indicating the characteristic energy scale which 
is lower for [110] than the other two, consistent with the 
different AF exchange in our model. While the shape of 
the susceptibility for [110] and [001] are very similar the 
latter direction is enhanced by a constant. This can be 
attributed to the temperature-independent orbital Van 
Vleck contribution along the z-axis, related to the or- 
bital operator r^. 

Lastly we consider the end member Ca2Ru04. The 
structure of Ca2Ru04 is orthorhombic (space group 
Pbca) [Q and is obtained from the high-symmetry tetrag- 
onal lA/mmm by rotating the RuOg-octahedra around 
their long axis (001) and tilting around the diagonal in- 
plane axis (110), with the condition that all planar oxy- 
gens move the same distance from the plane. As a result 
there is only one type of in-plane oxygen in this structure 
and all RuOg-octahedra are equivalent. 

Ca2Ru04 is an AF insulator and is relatively simple 
to describe with the LDA+U method jl^ . This method 
is based on spin-orbital unrestricted Hartree-Fock equa- 
tions (i.e. a static mean-field approximation) and gener- 
ally gives good results for long range ordered Mott insu- 
lators [^ . The Coulomb direct and exchange interaction 
parameters values used in our calculation were U=1.5 eV 
and J=0.7 eV. The groundstate of Ca2Ru04 was found 
to be an AF insulator with a sublattice moment of 1.51 
jj^B and a small energy gap, 0.30 eV, which compares well 
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Region 



I (0 < X < 0.2) 

II (0.2 < a; < 0.5) 

III {x -> 0.5) 
III (x = 2) 



orbital 



spin 



order 



(2,2) 
(3,1) 
(3,1) 



5 = 1 AF 
{yz,zx) S=\ FO, (AF) 
{yz,zx) S i (AFO) , (FM) 

S=0 



TABLE I. The three Regions of Ca2-^Sr^Ru04 with their 
orbital occupancy, the localized orbital and spin degrees of 
freedom and the order or (dominant correlations) (AF = spin 
antiferromagnetic, FM = spin ferromagnetic, AFO = anti- 
ferro-orbital, FO = ferro-orbital) . 



magnetic response. 
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with the experimental values 1.3 /is [0 and 0.2 eV |^ re- 
spectively. The band t2g is now fully spin-polarized (the 
reduced moment relative to S' = 1 is due to hybridization 
with oxygen states). Of the 4 t2g-electrons, 3 are in the 
majority-spin band and 1 is in the minority-spin band. In 
general, this implies an orbitally degenerate state. How- 
ever in Ca2Ru04, this degeneracy is lifted, since the xy- 
level lies considerably lower than the (xz, ?;z)-levels and 
is fully occupied in the insulating state. In the notation 
introduced above, we have integer (2, 2) occupancy of 
the orbitals in contrast to the (3,1) occupancy found for 
Regions III and II with x > 0.2. Using second order per- 
turbation theory for spin-orbit coupling, we find the mag- 
netic moment along the 6-axis in orthorhombic notations 
(or [110] direction in a tetragonal notation), in agree- 
ment with the experimental data |Q. The AF insulating 
phase extends to a; = 0.2 where by a first order structural 
change occurs to Region II (0.2 < x < 0.5) which exhibits 
short-range AF correlation and is metallic [^ . This first 
order boundary marks the transition from the (2, 2) to (3, 
1) configurations with increasing x. In summary, we have 
presented a consistent picture for the unusual phases of 
the isoelectronic alloy series Ca2-2:Sr2Ru04 based on the 
full multi-band electronic structure (see Table I). Start- 
ing from the good metal Sr2Ru04, we find the effect of 
Ca-substitution is to transfer electrons from the wider 
(xy)-band to the narrower (xz, yz) bands until at a crit- 
ical value of Xc = 0.5 there is integer occupancy of both 
subbands. The progressive rotation of the RuOg octa- 
hedra in this region leads to Mott localization of the 3 
electrons in the narrower (xz, yz) bands while the wider 
xy-band which is now half-filled, remains metallic. This 
partial localization of the 4d electrons can explain the 
puzzling observation of the coexistence of free S* = 1/2 
local moments and metallic behavior in Cai.5Sro.5Ru04. 
The ordered phases appearing in the Ca-rich Regions I 
(0 < X < 0.2) and II (0.2 < a; < 0.5), are based on local- 
ized (?/z, za;)-orbitals, while the xy-ovhital plays the role 
of a charge reservoir, yielding a completely filled (insu- 
lating)or a half-filled (metallic) band in the two regions, 
respectively. In the former case the ordered local spin 
has S* = 1, but in the latter case short-range correla- 
tions of a local S* = 1/2 spin combined with orbital order 
are realized, generating a pronounced anisotropy in the 
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